Three-parton contribution to pion form factor in kx factorization 
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We set up a framework for the study of the power-suppressed three-parton contribution to the 
pion electromagnetic form factor in the fcr factorization theorem. It is first shown that the gauge de- 
pendence proportional to parton transverse momenta from the two-parton Fock state and the gauge 
dependence associated with the three-parton Fock state cancel each other. After verifying the gauge 
invariance, we derive the three-parton-to-three-parton fcr-dependent hard kernel at leading order 
of the coupling constant, and find that it leads to about 5% correction to the pion electromagnetic 
form factor in the whole range of experimentally accessible momentum transfer squared. This sub- 
leading contribution is much smaller than the leading-order twist-2, next-to-leading-order twist-2 
' and leading-order two-parton twist-3 ones, which have been calculated in the literature. 
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I. INTRODUCTION 

O |. Aspects of the fcr factorization theorem [H-Q in perturbative QCD have been investigated intensively. One of the 
' ■ important issues is about the derivation of a fcy-dependent hard kernel at subleading level, which is defined as the 
difference between QCD diagrams and effective diagrams for transverse-momentum-dependent (TMD) hadron wave 
(— I functions. We have explained that partons in both sets of diagrams should remain off mass shell by in the kr 
I— factorization theorem |7|. The same statement has been made in the application of the fcr factorization theorem 
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PACS numbers; 12.38.Bx, 12.38. Cy, 12.39.St 



to inclusive processes such as prompt photon production [8[ . The ofF-shellness of partons may cause concern of the 
gauge invariance [1, However, we have shown that the gauge d ependence cancels between the above two sets of 
(yQ , diagrams, and a fcr-dependent hard kernel is gauge invariant [lil, [1^. Following this prescription, the next-to-leading- 
^\ ' order (NLO) correction to the pion transition (electromagnetic) form factor associated with the process ttj* — j{tt) 
on , has been calculated at leading twist, i.e., twist 2 [I3g Here we shall study the power-suppressed three-parton 
' contribution to the pion electromagnetic form factor in the fcr factorization theorem. The three-parton contribution 
' in the collinear factorization theorem [15^ to a simpler process, the p meson transition form factor, has been evaluated 
I recently fl6| . 

We shall first demonstrate the gauge invariance of the three-parton contribution to the pion electromagnetic form 
factor in the kx factorization theorem. There are two sources of gauge dependence for this power correction [l7j : 
the first source is proportional to parton transverse momenta from the two-parton Fock state. The corresponding 
hadronic matrix element is written as 



(0|g(z)n9„(7(0)|7r), (1) 

where z is the coordinate of the anti-quark field q, and F represents a combination of Gamma matrices. The second 
source is associated with the three-parton Fock state with an additional valence gluon. The corresponding matrix 
element is given by 

(0|g-(z)Fgr«A«(z')g(0)k), (2) 

with a color matrix T°, and the gluon field A° at the coordinate z'. The gauge dependences from the above two 
sources cancel each other, when Eqs. ([T]) and ^ are combined to form the gauge-invariant matrix element 

{0\q{z)TiD^{z')q{0)\Tr), (3) 
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^ Criticisms raised in ,10,1 have been responded in [Til . [T^ . 
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with the covariant derivative iDa = ida + qT°'A'^, . The cancellation of the gauge dependences is similar to that 
occurring in the collinear factorization theorem |17| . In the kr factorization theorem we just keep the transverse 
momentum dependence in denominators of particle propagators Hence, it is natural that the gauge dependence 
disappears at higher twists in the same way as in the collinear factorization theorem. 

Our formalism implies that contributions proportional to transverse momenta in numerators of hard kernels must 
be combined with contributions from three-parton Fock states in order to guarantee gauge invariance. Therefore, 
there is concern on the study of the pion transition form factor in [l^, where only the former was included. It has 
been pointed out explicitly that the leading-order (LO) hard kernel for the pion electromagnetic form factor becomes 
gauge-dependent, if one simply considers parton transverse momenta in numerators [loj . However, the contribution 
from the three-parton Fock state was still missing in , such that the false postulation on the gauge dependence of 
the kx factorization theorem was made. 

Since both the initial- and final-state pions involve higher- twist matrix elements like that in Eq. ([3]), the three- 
parton contribution to the pion electromagnetic form factor is suppressed at least by being the momentum 
transfer squared. After examining the gauge invariance, we calculate the fcr-dependent hard kernel for the three- 
parton-to-three-parton scattering in the Feynman gauge, and convolute it with the three-parton pion wave functions. 
We observe that the diagrams with a four-gluon vertex dominate this power correction. It will be shown that the 
three-parton contribution is only about 5% of the sum of those which have been investigated before, including the LO 
twist-2, NLO twist-2, and LO two-parton twist-3 ones fT3|. That is, the three-parton contribution is not crucial for 
accommodating experimental data of the pion electromagnetic form factor. At the same power of one should 
also take into account the scattering of two (four) partons into four (two) partons in principle. This piece has been 
analyzed in light-cone sum rules [T^, and found to be less important than other contributions. With this work, we 
conclude that the chirally enhanced two-parton twist-3 correction is the most important correction to the pion 
electromagnetic form factor. 

In Sec. II we verify the gauge invariance of the LO three-parton contribution to the pion electromagnetic form factor 
by combining the gauge-dependent hard kernels corresponding to Eqs. ([1]) and ©. The fcT-dependent hard kernel is 
then derived from the three-parton-to-three-parton scattering diagrams, and convoluted with the three-parton pion 
wave functions numerically in Sec. HI. Section IV is the conclusion. Detailed calculations of the gauge-dependent hard 
kernels corresponding to Eq. ([2]) are presented in Appendix A, and the expressions of the three-parton-to-three-parton 
hard kernels are collected in Appendix B. 



II. GAUGE INVARIANCE 



Consider the pion electromagnetic form factor involved in the process 7r(Pi)7* — )■ 7r(P2), whose LO diagrams are 
displayed in Fig. [TJ The momentum Pi of the initial-state pion and P2 of the final-state pion are parameterized as 

Pi = (P+,0,Ot)--^(1,0,Ot), P2 = (0,P2",Ot) = ^(0,1,Ot), (4) 

with — —q^, q ^ P2 ^ Pi being the virtual photon momentum. The gluon propagator of momentum / is written as 



'->'-). (^) 

in the covariant gauge, where the parameter A will be used to identify sources of gauge dependence. We assume that 
the anti-quarks in the initial- and final-state pions, represented by lower fermion line, carry the parton momenta 

fcl = (.TiF+,0,kiT), fc2 = (0,X2P2",k2T), (6) 

respectively, xi and X2 being the momentum fractions. It is understood that the components ki and A:^ have been 
dropped in hard kernels, and integrated out of the TMD pion wave functions. 
We employ the Fierz identity 



lijhk 



to factorize the fermion fiow, where / denotes the 4x4 identity matrix. The structure 7q75 in the above identity 
contributes at twist 2 and higher twists, and 75 and CTq,^75 contribute at twist 3 and higher twists at two-parton level. 
The matrix F in Eqs. ([T]) and ([2]) can pick up one of the above structures, among which we focus on the first one 7^75 
as an example below. We also insert the identity 

hjhk = j^kjlik + 2{T^')ij{T'')ik, (8) 
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(a) (6) 
FIG. 1: LO diagrams for the pion electromagnetic form factor, where the symbol x represents the virtual photon vertex. 

to factorize the color flow, where A^c = 3 is the number of colors, / denotes the 3x3 identity matrix, and T"^ is a color 
matrix. The first (second) term in Eq. (|5]) will be associated with a color-singlet (color-octet) state of the valence 
quark and anti-quark. 

The first source of gauge dependence is extracted from the diagrams in Fig. [U where the quark and anti-quark pair 
forms a color-singlet state. Combining the decompositions in Eqs. ([7]) and ([HI), we sandwich Fig.[T]with the structures 

^7.75, ^757/^, (9) 

from the initial and final states, respectively, where the subscripts a and /3 can take arbitrary components. The LO 
hard kernel from Fig. [IJa) contains the gauge-dependent piece 

rraX^ - 2. ^?-[7'^757/37m ( /^l " hh^lalb] {kl - fc2)a(fcl - fcz)^. , . 

'^^^16N, {Pi-k^y{k,-k^Y {k,-k^Y ' ^'"^ 

with Cp being a color factor. The above expression diminishes with the substitution ki = xiPi, fca = X2P2, 7q — 7^ 
(proportional to pi), and — 7+ (proportional to ^2) into the numerator, implying that Eq. ([T0|) does not contribute 
at leading power in the kx factorization theorem. 

To obtain the gauge-dependent hard kernel from Eq. ([TU|) at l/Q^, we insert the identity (fci — k2)i, — (Pi — fc2)iy — 
(Pi — ki)^. It can be shown that the contribution from the (Pi — /c2)i/ term is cancelled by the corresponding one 
in Fig. [ijb). The second term, with Pi — ki being the momentum of the incoming valence quark, corresponds to 
the matrix element with the derivative of the quark field in the initial-state pion. This term can be picked up by 
differentiating Eq. ((TU)) with respect to ki . Once Eq. ([T]) for the initial-state pion is identified, we pick up the k2a term 
in (kl — ^2)0- via differentiation, which corresponds to the derivative of the valence anti-quark field in the final-state 
pion. Note that denominators of particle propagators, depending on ki and k2, will be differentiated too. However, 
their differentiation gives rise to even higher-twist matrix elements, and can be neglected. We then extract 

HrAki,k2} - g^^^g^^^ - ^eg X ^^^^ ^^^^^^^ , (11) 

associated with (0|g(2;)757"i(9Q(7(0)|7r(Pi)) and the similar matrix element for the final-state pion. 
The LO hard kernel from Fig. [Hb) contains 



rrbX^ - 2. tr[j''-f5-fi3Yi^2- Mht^lalb] (kl - fc2)a(fcl - k2)^ 

16iV, (p^-kif{ki-k2f {kl-k2Y ' ^ ' 



whose differentiation with respect to fcic and /c2/3 leads to 

U^>^(h h\ - .■-.-2^ Cf tr-[7"757ff7^(/P2-a;i /Pi)7^7c.75] , „x 

Hj,j,{ki,k2) - teg A—— — j—rjTT T^a ■ (l-^i 

16Nc {P2- kiY[ki- k2r 

This hard kernel corresponds to the matrix element of the initial-state pion with the derivative of the anti-quark field 
g(z). Equations (jlip and ()13p represent the gauge dependence in LO two-parton-to-two-parton scattering at power 
of l/Q^, which was also observed in [Toj . 

The second source of gauge dependence arises from Fig.[2l where three partons appear in the initial state as indicated 
by Eq. All possible attachments of the additional valence gluon to the lines other than the valence quark and 
anti-quark in the initial state are labeled by letters A, B, ■ ■ ■ . In principle, the diagrams with the attachments to 
the valence quark and anti-quark in the initial state should be included in order to respect U{1) gauge symmetry for 
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FIG. 2: Diagrams with three partons from the initial state, where letters A, B, ■ ■ ■ denote the attachments of the additional 
valence gluon. 



the electromagnetic interaction^. These diagrams contribute to higher Gegenbauer terms in the two-parton twist-3 
pion distribution amplitudes. Equations of motion can then be constructed to relate the coefficients of the higher 
Gegenbauer terms in the two-parton twist-3 and three-parton twist-3 pion distribution amplitudes pol |. Hence, 
one should pay attention to the consistency between the models for these two sets of distribution amplitudes in a 
numerical analysis. We shall adopt the non-asymptotic models for both sets of distribution amplitudes in Sec. Ill, 
when estimating the importance of the three-parton contribution relative to other two-parton contributions in the 
pion electromagnetic form factor. 

According to Eq. ([8|), we absorb a color matrix T'^ and the coupling constant g associated with an attachment 
into the matrix element for the initial-state pion, and another T'^ goes into the evaluation of gauge-dependent hard 
kernels. For example, the color factor corresponding to the attachment A of the valence gluon to the virtual quark 
line is given by 

where the color matrix T° (T^) comes from the valence (hard) gluon vertex. After summing over c, the tensor S"''^ 
sets c a in the matrix element for the initial-state pion, leading to Eq. Including the coefficient 2 in Eq. we 
adopt the structure 7^75/2 for the initial state in the calculation of Fig. [21 whose details can be found in Appendix 
A. The results are collected as follows: 



327V2 {Pi-k2nh-k2) 



rX _ . 2x 1 ^^[7'^757/37m(/'i - 2^2 Ah^lalb] 



""^^ ^ '32 ip,-k2nk,-k2f ' 



H^T = 0, (17) 



Hut - ~^eg (P, - - ^^(fc, - fc,)4 ' (^8) 



r\ _ ■ 2. ir[7^757/j7"(/p2 - a;i a - yi /F'i)7A*7a75] 



""-^ = ""^ (P2 - fci - hnk, - k2Y ' ^''^ 



1 [7^757^7" (/P2 -X2 A- yi /Pi)7M7a75] 
327V| (P2_A:2-/i)2(fci-fc2)4 



. 2x 1 ^?-[7^757/37"(/2 -xi A- yi /Pi)7^7a75] , . 



327V2 ^p^_ki-hY{ki-k2Y 
1 ^?-[7^757/37"(/2 -xi Pi- yi A)7A»7a75] 



1 fr[7"757^7^(/p2 - xi /Pi - yi /Pi)7^7q75 



^ 32 {P2~ki~hy{ki-k2r 

A— — , , , rri , (21) 



32 (P2-fci-/i)2(fci+/i 



rA _ ■ 2\ 1 ^''[7"75707'^(/2 - a;i /Pi - yi /Pi)7^7a75] 



- ^32N! {P2-ki-h)Hki+h-k2r ' ^ ' 



^ We thank V. Braun for pointing out this U{1) gauge symmetry. 
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FIG. 3: One of the three-parton-to-three-parton diagrams, where the two valence ghions scatter via a three-gluon vertex. 



with the gluon momentum fraction yi = / . 
Summing the above expressions, we arrive at 



H 

Y^H^r^H^gT + Hfr, (23) 

i=A 



with 



- i/^^(fci, fc2), if„T - H^riki + WM)- (24) 



The contributions from the attachments A and B are added into H^^ with the desired color factor Cp- The con- 
tribution from the attachment D and the first term in the attachment F cancels each other. The second term from 
the attachment F and the first term from the attachment G are combined into the expression with the color factor 
Cp, which then cancels the contribution from the attachment E. The second term of G and the contribution from 
the attachment H are added into FI^t- Note that H^^ does not depend on the valence gluon momentum li, which 
can then be integrated out of the matrix element, giving {Q\q{z)^]^^"-gT°-A'^{Q)q{Q)\-K{Pi)). The hard kernel -ffgj-; 
depending on the combination fci + /i, corresponds to the matrix element {Q\q{z)"f^"f°' gT°- A'^{z)q{Q)\'K{Pi)) . Because 
of the symmetry under the exchange of the initial- and final-state kinematic variables, the gauge-dependent hard 
kernels with three partons from the final state are written as 

Ht) = H^Uki, k2 + h), H^r'a - Hf^ikiM), (25) 

where I2 is the momentum carried by the outgoing valence gluon. 

At last, we extract the gauge dependence from the three-parton-to-three-parton diagrams. In this case the indices 
a and /?, associated with the initial and final valence gluons, respectively, must be carried by gluon vertices, instead 
of by parton momenta. Focusing on the gauge-dependent piece, we can apply the Ward identity to hard gluons. It is 
easy to find that Fig. [3l where the two valence gluons scatter via a three-gluon vertex, does not contribute: if the gauge 
dependence arises from the lower hard gluon, the results, being proportional to the parton momenta after applying 
the Ward identity, should be dropped. If the gauge dependence arises from the upper hard gluon, the Ward identity 
diminishes the amplitude for a similar reason. It is also easy to see that Fig. IHJa) with a four-gluon vertex does not 
contribute to a gauge-dependent hard kernel. If both the valence gluons attach to the quark line, the gauge-dependent 
contribution vanishes because of the Ward identity applied to the anti-quark line. If both the valence gluons attach 
to the anti-quark line, the gauge-dependent contribution vanishes too. 

The other diagrams are classified into several sets, which are formed by all possible attachments of a hard gluon as 
displayed in Figs. |4l[a)-|4l[l). Cancellation occurs in each set of diagrams after applying the Ward identity. Neglecting 
those pieces proportional to the parton momenta, we have 

(a) + (6) = (c) + {d) - 0, 

(e) + (/) + Cg) + {h) = 0. (26) 
A finite gauge-dependent contribution comes only from Figs. |4l[i)-|4ljl), given by 

(k) + (/) = H'g^ = Hl^^^ik, +h,k2). (27) 



(a) 



(b) 



id) 




(/) 



(9) 






FIG. 4: Three-parton-to-three-parton diagrams, where the valence gluon from the initial state is shown, and the possible 
attachments of the valence gluon from the final state are represented by dots. 



Due to the same expressions of Hf^ and ff^y, their corresponding matrix elements are combined into 

(0|g(z)75»7"^a(0)g(0)k(Pi)) = 0, (28) 

with the equation of motion for the quark field, /D(0)q(0) = 0. That is, the gauge invariance holds, when the 
contributions from Figs. [T] and [5] are combined. The combination of the matrix elements for iJ^^ and H^g also 

vanishes according to Eq. ([28| . A similar reasoning applies to the gauge-dependent hard kernels and and 
to H^gT fl'iid H^g'. the combination of their matrix elements vanishes due to the equation of the motion for the quark 
field in the final-state pion. This observation completes the proof of the gauge invariance of the LO three-parton 
contribution to the pion electromagnetic form factor at power of 1/ in the fc-p factorization theorem. The extension 
of the proof to all orders can follow the steps outlined in [13,] . Note that our proof applies to the collinear factorization 
theorem too: simply neglecting transverse momenta in denominators, one can show the gauge invariance of the 
three-parton contribution to the pion electromagnetic form factor in the collinear factorization theorem. 



III. THREE-PARTON CONTRIBUTION 



In this section we calculate the three-parton contribution to the pion electromagnetic form factor. Start with the 
gauge-invariant twist- 3 matrix element 

(0|q-(z)a+„,75»^a(^')'Z(0)|vr(Pi)), (29) 

where the subscript a is associated with the vertex the valence gluon attaches to. The power behavior will not 
be changed, and the gauge invariance will not be broken by inserting another covariant derivative 13"^. We then 
exchange -D+ and Da, take the difference of D^Da and DaD^, and apply the identity [D^,Da] = —igG^^- It is 
then equivalent to employ the following alternative matrix element 21], which defines the three-parton twist-3 pion 
wave function T{z, z'), 



(0|g(z)a+„,75ffG+Jz')9(0)k(Pi)) = iUmo{P+)'gla'Tiz, z'), 



(30) 
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with the chiral scale uiq = rn^/iniu + md), rriT^, rriu, and being the pion, u quark and d quark masses, respectively. 
The operators with other spin structures contribute at higher twists: for example, the operator ^^"y^Ga/j gives a 
three-parton twist-4 contribution, and ^^Gaf) does not contribute [13, US 113 ■ With Eq. (pO| . one may verify the 
gauge invariance of a hard kernel in the fc^ factorization theorem by demonstrating the cancellation between the 



gauge-dependent contributions from the operators and daA'^ 21 [ 

Below we derive the hard kernels from the three-parton-to-three-parton diagrams corresponding to Eq. pop in the 
Feynman gauge (A = 0). Choosing this gauge, the operator 9qA+ does not contribute, so only is relevant. The 

three momenta Pi — fci — ^i, fci, and are assigned to the initial-state quark, antiquark, and gluon, respectively, and 
P2 — k2 — l2,k2, and I2 to the final-state quark, anti-quark, and gluon, respectively. We have the structures for the 
initial- and final-state pions 

4 1+ Ayi 

-l5^cr^'^' j:^hi)f7Tmo{P^fgJf^, = /p27j/7rWo, (31) 

where the gluon momentum fraction j/2 is defined by 1/2 = I2 1-^2 ^ t^*^ gamma matrix 7-^ involves only transverse 
components. 

There are totally 196 diagrams for the three-parton-to-three-parton scattering, which can be divided into four 
categories'^. Category A contains 20 quark-gluon configurations, in which neither of the valence gluons attaches to 
the hard gluon line. Each configuration allows 6 different attachments for the photon line. We further divide this 
category into two groups as shown in Fig. [51 where both valence gluons attach to the same quark or anti-quark 
line, and in Fig. [71 where one valence gluon attaches to the quark line and another to the anti-quark line. Only the 
diagrams giving nonvanishing contributions are displayed. As observed in Appendix B, the amplitudes with both the 
valence gluons attaching to the quark line are power-suppressed. Category B contains 8 quark-gluon configurations, 
in which one of the valence gluons attaches to the hard gluon line. Each configuration allows 5 different attachments 
for the photon line, among which those with nonvanishing contributions are displayed in Fig. [51 Category C contains 
4 quark-gluon configurations, where both the valence gluons are connected to the hard gluon. Each configuration 
allows 4 different attachments for the photon line, among which those with nonvanishing contributions are displayed 
in Fig. [31 Category D contains 4 quark-gluon configurations, where the two valence gluons scatter via a three-gluon 
vertex as shown in Fig. [31 Each configuration allows 5 different attachments for the photon line. Since this category of 
diagrams does not contribute, we shall not discuss them further. Besides, when a valence gluon attaches to a valence 
quark, the diagram should be regarded as being from an effective two-parton Fock state, and will not be calculated. 

We extract the hard kernels proportional to the final-state momentum P2^. The hard kernels proportional to Pi^ 
can be obtained by exchanging the kinetic variables of the initial- and final-state pions. Adopting the electric charge 
e, instead of the quark charge e„ or ea, we have taken into account the diagrams with the virtual photon attaching 
to the anti-quark line. Figure [SJa) with a four-gluon vertex gives the dominant three-parton contribution 

16yi?/2 8(N^ - 1) (Pi - k2- hriki - k2Y{ki - k2 + h - I2Y 

^ _. 2 />o tr[/P2^, /P2 /Pi] 

^ lQyiy2 - 1 (Pi - P2)Hk, - k2 + h - h^iki - fca)^ ' ^ ' 

To arrive at the second line, we have made an approximation according to the power counting 3> xQ^, yQ^ ^ fc|, 
[l^ . under which the TMD term in the following denominator is neglected, 

(A2+/2) _ {k2+l2)-l+ _ P2I+ _ A ^33^ 



{P^-k2-l2f ~2P+{k2+l2)- -"^P^Pi {Pl-P2f 

The expressions for other three-parton-to-three-parton diagrams are collected in Appendix B. 

Since the Sudakov factor for exclusive QCD processes was derived in the space of impact parameters 0, [H| , we 
Fourier transform Eq. p2|) . The fcr factorization formula for the pion electromagnetic pion form factor from Fig.[9l^a) 



^ We thank the Referee for suggesting this classification of diagrams. 
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is then written as 



dxi I — ^ / dx2 
Vi Jo 



y2 



x<^{xi,yi)^{x2,y2)K [y^{xi+yi){x2 + 2/2)Q K{y/^^Q), 



(34) 



where the three-parton pion distribution amphtude $(a;i,yi) corresponds to T(z,z') in Eq. pop in the space 
momentum fractions. The functions K, arising from the Fourier transformation of the TMD denominators (fci — /c2 
h — hY and {ki — k2Y in Eq. ([5^. are defined by 



of 



Kit) = 



l/A 



bdbKoitb) cxp[-s{P+,b)], 



(35) 



in which Kq is the modified Bessel function, and the exphcit expression of the Sudakov exponent s{P^ , b) is referred 
to [1^ . We have kept only the most effective piece of the Sudakov evohition in the smah x region, that results from 
the gluon exchanges between the energetic valence quark and the Wilson line associated with it. Because the Sudakov 
factor exp[— s(P{'', 6)] diminishes at & = 1/A, with the QCD scale A ss 0.3 GeV, the upper bound of the integration 
variable has been set to 1/ A in Eq. (|35|) . For order-of- magnitude estimate and for demonstrating the smallness of the 
three-parton contribution, we do not consider the renormalization-group evolution from the low scale, at which $ is 
defined, to the scale of the hard kernel. The coupling constant is also assumed to be a constant — 0.5. 

Fourier transforming the hard kernels in Appendix B into the impact-parameter space, we construct the corre- 
sponding kT factorization formulas similar to Eq. (I34[) . We then add all the contributions to the pion electromagnetic 
form factor, and employ the model of the three-parton twist-3 pion distribution amplitude (20, 22] 



$(xi,yi) = 360773a;i(l ~ xi ~ yi)yl 



i + f(7yi-3) 



(36) 



for a numerical analysis, with the parameters 773 = 0.015 and = —3, and 1 — xi — yi being the momentum fraction 
of the valence quark. The total three-parton contribution to Q^F{Q^), F{Q^) being the pion electromagnetic form 
fact, is displayed in Fig. [S] The curve exhibits a decrease in (though not obvious in the figure) compared to the 
LO and NLO twist-2 contributions, indicating that this contribution is power-suppressed. It is only about 5% of the 
sum of those evaluated in [3 HE], including the LO twist-2, NLO twist-2, and LO two-parton twist-3 pieces. That 
is, the three-parton contribution is not crucial for accommodating the experimental data of Q^F{Q^) [2^. [27j in the 
whole accessible range of up to 10 GeV^. The only important subleading contribution to the pion electromagnetic 
form factor that have been investigated so far comes from the chirally enhanced two-parton twist-3 one. 



IV. CONCLUSION 



In this paper we have applied the fcr factorization theorem to the study of the power-suppressed three-parton 
contribution to the pion electromagnetic form factor. It was demonstrated that the gauge invariance of the kr- 
dependent hard kernel holds for this power correction: the gauge dependence proportional to parton transverse 
momenta from the two-parton Fock state and the gauge dependence associated with the three-parton Fock state 
cancel each other. We have calculated the three-parton-to-three-parton hard kernel at LO, and found that the three- 
parton contribution is about 5% of the sum of the LO twist-2, NLO twist-2, and LO two-parton twist-3 ones in the 
whole range of experimentally accessible . Our analysis shows that the power expansion for this exclusive process 
might be reliable in the kx factorization theorem. At the same power of l/Q^, the two-parton twist-4 contribution 
should be taken into account, which has been studied in the framework of light-cone sum rules We shall calculate 
this correction in the kr factorization theorem in the future, which involves a twist-2 distribution amplitude from 
one side and a two-parton twist-4 distribution amplitude from the other side. We shall also extend our framework to 
exclusive B meson decays, for which three-parton contributions have been analyzed in light-cone sum rules, [28j . in 
the QCD (coUinear) factorization ^] and in the soft-coUinear effective theory |30 |. 

We thank Z.T. Wei for his suggestion, which initiates this project. We also thank V. Braun, S. Brodsky, J. P. Ma 
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FIG. 5: dependence of the three-parton contribution to the pion electromagnetic form factor. The LO twist-2, NLO twist-2, 
and LO two-parton twist-3 contributions are quoted from Fig. 8 in 14], which were derived using the non-asymptotic models 
for the pion distribution amplitudes. 



Appendix A: GAUGE DEPENDENCE 

In this appendix we present the detailed derivation of the gauge-dependent hard kernels from Fig. [5] The attachment 
A contains, with the color factor in Eq. (|T4l) . the gauge-dependent piece 



Inserting the identity — {/Pi— /i— 7^2) ^ (/Pi^ /i— M) i^^to Eq. (jAip . we obtain 

tr[{M- h)Toli3lt,{^i- h2)T{IP\- A- /:2)(/Pi- A- /l)7a75] 1 



(Pi-k2f{Pi-h-k2Y(ki-k2Y 



(A2) 
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The second term, proportional to the momentum Pi — ki — li of the incoming valence quark, should be dropped, since 
there is the valence gluon A" from the initial state already. Taking the derivative of the first term with respect to 
^2/3, and then substituting fci = 2^1 -Pi and ^2 = 2^2 ^2 into the numerator, we have Eq. (|15p . 

The diagram with the attachment B of the valence gluon to the hard gluon line produces the gauge-dependent hard 
kernel 

-"B — — tr{i 1 1 )Lg 



[Pi - 



(fcl - fc2)2 



[ki~k-,+hY{ki-k2+lir' , Aki-k2)'(ki-k2f 

V A- 



(fcl - + hY 



(ki - k2Y 



(fcl 



-A 



2 (fcl - fc2)*(fcl - fc2)^' (fcl - fc2 + /l)"(fcl - fc2 + hY 



(fcl - kY^ 



■CI-K2)- iki-k2 + hY 

with the triple-gluon vertex, 

pdfca ^ fdba[g^^(^2h + fcl " k2)s + 9us{2k2 ~ 2fci - hY + 5fc(fcl " fc2 " ^l).] 

jdba ^jgjj^g antisymmetric tensor. Using the identity 



^dbc ]-,gij-^g a symmetric tensor, the above amplitude becomes 
2 1 ir[75/757;37;,(/Pi- /:2hi^'l"l5] 



dbc £dba 



0, r^f 



(A3) 

(A4) 
(A5) 



'"^ 32 (Fi - k2Y 

X[ga„{2li + fcl - k2)s +gu5{2k2 - 2fci - li)a + gSa{kl - fc2 

9 



(ki^k2+hYiki-k2 + hY' . , (fcl - fc2)''(fci - 



+ A- 



fc2)'" 



(fcl-fc2)2 (fcl-fc2+/l)4 (fcl-fc2)4 

(fcl - fc2)*(fcl - fc2)'" (fcl - fc2 + ;i)"(fcl - fc2 + hY' 



(fcl - fc-2 + hY 



(fcl - fc2)^ 



(fcl - fc-2 +Zl)^ 



(A6) 



The term leads to 



2 1 \2^^[(A- /fc2)757/37M(A- A2+ /l)7a75] 



ieq — A 
^ 32 



(Pl-fc"2)2(fcl-fc2)4(fcl 



hY 



X (fcl - fc2) • (2/1 + fcl - fc2)(fcl - fc2 + hY + (fcl - fc2) • (fcl - fc2 + Zl)(2fc2 - 2fci - ^i)' 
+ (fcl - fc2 - h) ■ (fcl - fc2 + Zl)(fcl - fc2)" . 



(A7) 



Inserting the identity /;i— /:2+ /i = (/Pi— /:2) — (/Pi^ Ai^ A); if is easy to see that the first term of the identity 
gives an expression which is cancelled by the corresponding one from the attachment G. The second term, being 
proportional to the momentum Pi — fci — h of the incoming valence quark, should be neglected. 
We then consider the A terms. The first A term in Eq. (|A6[) gives 

1 tr[{2 /i+ /ni- /;2)757/37/x(A- h+ /i)7a75] 



" 32-^ (Pl-fc2)^(fcl-fc2)^(fcl-fc2 + /l)4 " + 

■ 2 1 X ^''[(A- /fc2+ /l)757/37p(/Pl- /fc2)(A- /fc2+ /l)7a75] „, 

32" (Pi-fc2)^(fci-fc2)^(fci-fc2+^l)^ (''^ - - 

■ 2 1 , ^r[7"757^7^(/Pi- h+ hhal5] n u i 1 \ (i, i. ; ^ 

-''^ 32^ (Pl-fc2)^(fcl-fc-2)^(fcl-fc2+^l)^ '^'^ ■ 



(A8) 



which are all negligible for the same reason as for the A'^ term. Hence, receives a contribution only from the 
second A term. 



^ '""^ 32^ (Pi - fc2)2(fci - fc2)4(fci - fc2 + hY 

2 1 .HiM- MhblPlt^i^l- /fc2)7"7a75] 



(fcl - fca) • (2^1 + fcl - fc2) 



32^ 



(Pl-fc2)2(fcl -fc2)4 



(A9) 
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FIG. 6: Three-parton-to-three-parton diagrams in Category A, where both valence gluons attach to the same quark or anti- 
quark hne. 



To arrive at the second line, the higher-power term = —l^x been added, so we have (fci — fc2) • (2/i + ki — k2) + lf = 
(fci — k2 + h)'^ ■ The differentiation of the above expression with respect to fc2/3 leads to Eq. (IT51) . The diagrams with 
other attachments in Fig. [5] can be calculated in a similar way, so the detail will not be presented here. 



Appendix B: HARD KERNELS 



In this appendix we collect the expressions of the three-parton-to-three-parton hard kernels for the pion electro- 
magnetic form factor in the Feynman gauge. Start with Category A defined in Sec. III. When the attachments of the 
two valence gluons are arranged in the way that the hard gluon vertices sandwich the spin structures associated with 
the pions, the contribution diminishes because of Pill^lv = or /27J71/ = 0. The nonvanishing amplitudes 
come from Figs. [H] and [71 which are written as 



„ ^ 2 ^K/P27p /P2 A A /Pi) , . 

N^{N^ - 1) (Pi - P2)2(Pl - fc2 - h^iPl ~k,-h~ P2)2(fci - ^2)2 ' ^ > 

rr ^ ■ 2 Nl + l ^r(/p27M A /Pi h A) 

" N^{N^-i){p,-p,)2(k,~k2 + h-i2r{ki-h)Hk2-h)^' ^ ^ 

„ ^ ■ 2 ^c+l tri^2 A /P27m /P2 /Pi) 

N^{N^ - 1) (P2 - A;2 - ^2 - Pl)2(P2 - fci - ;i)2(Pi -k^-W- P2)2(fci - k2Y ' ^ ' 

rr ^ ■ 2 1 tr{/P2 /l /p27M ^2 A) fj... 

Ar2(7V2_i)(p^_A:i-/i-P2)2(fci-fc2 + ;i)2(/i-fc2)2(P2-fci-;i)2' ^ ^ 

rr ^ _,,„2^ H^2l, A A /P2 A) 

" ^ N^AN^ - 1) (Pi - P2)2(fci - ^2 + h)Hh - fc2)2(Pl - fci - /l - P2)2 ' ^ ^ 

„ ^ . 2 1 tr(A7M A A A A) .pj.x 

'"'^ 7V2(7V2 _ 1) (p^ _ P2)2(fc^ _ ^2)2(fci - Z2)(Pl - fc2 - h)^ ' ^ ' 

It is observed that the results from Figs.[6ja) and[6jc) are suppressed by a power of l/Q compared to Fig.[6][b). That 
is, when both the valence gluons attach to the quark line, the contribution is power-suppressed. The hard kernels 
from Fig. [7] are of the same power as Eq. (|B2p . 
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FIG. 7: Three-parton-to-three-parton diagrams in Category A, where one valence gluon attaches to the quark hne and another 
to the anti-quark hne. 
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FIG. 8: Three-parton-to-three-parton diagrams in Category B, where one valence gluon attaches to the hard gluon line. 



Figure [5] from Category B gives the hard kernels 

1 tr[^2{M- A) ^2-/^ A /Pi] 



Hsa 

Hsb 

Hse 



2 ^ 7V2 _ 1 (p^ _ fci - -k2+ hYiki - k^YiPl -ki-h- P2f ' 



1 



tr[^2l^. A) A /Pi] 



2^-^ iV2 - 1 (Fi - P2)2(fci - + /l)2(fci - fc2)2(Fl - fci - - P2)^ 

i 2 tri^2l^ ^2 /Pi) 



2 (Pi-P2)2(fci-A;2 + ?i)2(fci-fc2)2' 



tK/P27M /P2 /Pl(/:2+2 A) A] 



2 ^ iV2 - 1 (Pi - P2)2(fci - ^2 + Zl)2(?i - k2)Hk, -k2+h- h? ' 

i_ 2_^ ^K/P27M /P2 /Pi A(A +2 A)] 

'2^^ 7V2 _ 1 (p^ _ P2)2(fc^ - + /l - Z2)^(fcl - /2)2(fcl -k2~ I2? ' 

i_ 2_}: ^K/P27p /P2 /Pi(^- h) A] 

2''^ 7V2 _ 1 (p, _ P2)2(fc^ _ fc2)2(fci - fc2 - /2)2(Pl - ^2 - /2)2 ' 



(B7) 

(B8) 
(B9) 
(BIO) 
(Bll) 
(B12) 
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FIG. 9: Three-parton-to-three-parton diagrams in Category C, where both valence gluons attach to the hard gluon hne. 



Figure [H] from Category C contributes the hard kernels 



9b 



-leg 



- 1 (Pi - P2)Hki -k2+h- h)^(h - fc2)2' 

, (fci-fc2-2/2 + ^i)-(fci-fc2-/i)tr(/p27^ /Pa A 



1 (-Pi - P2Y{kl -k2+ hYiki - fc2)2(fcl - fc2 + ^1 - k? ' 

Nl (fcl - fc2 + I2) ■ (fci - fc2 + 2h - Z2)tr(/P27M jP2 /Pi) 



1(^1 



P2?{ki -k2 + h- hriki -k2- kriki - k2f 



(B13) 
(B14) 
(B15) 



[1 
[2: 

[3; 

[4] 

[6: 
[7: 

[lo; 
[11 

[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22' 
[23^ 
[24 
[25| 

[26 
[27 
[28' 
[29 
[30 



S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. B 242, 97 (1990); Nucl. Phys. B366, 135 (1991). 
J.C. Collins and R.K. Ellis, Nucl. Phys. B360, 3 (1991). 

E. M. Levin, M.G. Ryskin, Yu.M. Shabelskii, and A.G. Shuvaev, Sov. J. Nucl. Phys. 53, 657 (1991). 
J. Botts and G. Sterman, Nucl. Phys. B225, 62 (1989). 

H-n. Li and G. Sterman, Nucl. Phys. B381, 129 (1992). 

T. Huang and Q.X. Shen, Z. Phys. C 50, 139 (1991); J.P. Ralston and B. Pire, Phys. Rev. Lett. 65, 2343 (1990). 

M. Nagashima and H-n. Li, Phys. Rev. D 67, 034001 (2003)^ 

A.V. Lipatov, M.A. Malyshev, and N.P. Zotov, arXiv: 1 102^1341 [hep-ph]. 

F. Feng, J.P. Ma and Q. Wang, Phys. Lett. B 674, 176 (2009). 

F. Feng, J.P. Ma and Q. Wang. arXiv :1103.1706 [hep-ph]. 
H-n. Li and S. Mishima, Phys. Lett. B 674, 182 (2009). 
H-n. Li and S. Mishima, Phys. Rev. D 80, 074024 (2009). 
S. Nandi and H-n. Li, Phys. Rev. D 76, 034008 (2007). 

H-n. Li, Y.M. Wang, Y.L. Shen, and H. Zou, Phys. Rev. D 83, 054029 (2011). 

G. P. Lepage and S.J. Brodsky, Phys. Rev. D 22, 2157 (1980). 

I.V. Anikin, D.Yu. Ivanov, B. Pire, L. Szymanowski, and S. Wallon, Nucl. Phys. B828, 1 (2010). 
J.W. Qiu, Phys. Rev. D 42, 30 (1990); J. Qiu and G. Sterman, Nucl. Phys. B378, 52 (1992). 

F. G. Cao, T. Huang, and B.Q. Ma, Phys. Rev. D 53, 6582 (1996). 

J. Bijnens and A. Khodjamirian, Eur. Phys. J. C 26, 67 (2002); S.S. Agaev, Phys. Rev. D 72, 074020 (2005). 
V.M. Braun and LE. Filyanov, Z. Phys. C 48, 239 (1990). 

H. Eguchi, Y. Koike, and K. Tanaka, Nucl. Phys. B763, 198 (2007). 
P. Ball, J. High Energy Phys. 01, 010 (1999). 

P. Ball, V.M. Braun, and A. Lenz, J. High Energy Phys. 05, 004 (2006). 

Y.Y. Keum, H-n. Li, and A.I. Sanda, Phys Lett. B 504, 6 (2001); Phys. Rev. D 63, 054008 (2001). 
Z.T. Wei and M.Z. Yang, Phys. Rev. D 67, 094013 (2003); J.W. Chen, H. Kohyama, K. Ohnishi, U. 
Phys. Lett. B 693, 102 (2010). 

G. M. Huber et al. [Jefferson Lab Collaboration], Phys. Rev. C 78, 045203 (2008). 
C.J. Bebek et al, Phys. Rev. D 17, 1693 (1978). 

A. Khodjamirian, Nucl. Phys. B605, 558 (2001). 

K.C. Yang, Phys. Rev. D 69, 054025 (2004); T.W. Yeh, larXiv:0802.1855l [hep-ph]. 

CM. Arnesen, I.Z. Rothstein, and I.W. Stewart, Phys. Lett. B 647, 405 (2007); Erratum-ibid. B 653, 450 (2007). 



Raha, and Y.L. Shen, 



